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SUMMARY 

Numerical solutions are presented for the laminar and turbulent boundary -layer 
equations for incompressible flows with separation and reattachment. The separation 
singularity is avoided by using an inverse technique in which the displacement thickness 
is prescribed and the pressure is deduced from the resulting solution. The turbulent 
results appear qualitatively correct despite the use of a two-layer eddy-viscosity model 
which is generally assumed appropriate only for mild-pressure-gradient flows. A new 
viscous-inviscid interaction technique is presented in which the inviscid flow is solved 
inversely by prescribing the pressure from the boundary -layer solution and deducing the 
new displacement thickness from the solution of a Cauchy integral. Calculations are 
presented using this interaction procedure for a laminar flow in which separation and 
reattachment occur on a solid surface. 

INTRODUCTION 

The development of theoretical prediction techniques for flows involving boundary- 
layer separation is of fundamental importance since separation is a common occurrence 
on most aerodynamic surfaces. Significant progress has been made in recent years for 
laminar supersonic flows which contain large viscous-inviscid interactions, such as those 
which occur at the point of incidence of a shock wave on a flat plate or in the vicinity of .a 
compression corner. Werle and Vatsa (ref. 1) and Dwoyer (ref. 2), as well as others, 
have demonstrated that the boundary-layer equations including interaction with the invis- 
cid flow provide an accurate model which gives results that agree with experiment and 
solutions of the Navier- Stokes equations (ref. 3). In addition these techniques are cur- 
rently being extended to turbulent flows (ref. 4). 

In supersonic flows the interaction between the viscous and inviscid flow can be 
computed locally since the inviscid flow is hyperbolic. As a result the singularity which 
occurs in solutions of the boundary-layer equations when the pressure is prescribed is 
removed by allowir^ the boundary layer to modify the pressure to give a regular solution 
at the separation point. This technique cannot be used for subsonic flows since the invis- 


125 


cid flow is elliptic and hence, the pressure at any given point depends on the entire dis- 
placement body distribution. 

Alternate procedures which are appropriate for subsonic flows and do not rely bn 
interaction to eliminate the separation singularity are the inverse boundary-layer solu- 
tion techniques which have been recently developed. Catherall and Mangier (ref. 5), and 
later Carter (ref. 6), have demonstrated that regular solutions can be obtained if the dis- 
placement thickness is prescribed and the pressure is deduced from the resulting solu- 
tion. Similarly, regular solutions at separation can be obtained by prescribing the skin 
friction as shown by Kuhn and Nielsen (ref. 7), Klineberg and Steger (ref. 8), and Carter, 
(ref. 6). Another inverse procedure, which was developed earlier by Klineberg and 
Steger (ref. 9) and later used by Tai (ref. 10), is to use the transverse component of 
velocity at the boundary-layer edge as the prescribed condition. 

In contrast with these inverse techniques, Briley and McDonald (ref., 11) have made 
calculations for subsonic flow using a direct procedure in which the unsteady boundary- 
layer equations are repeatedly solved until a steady-state solution is obtained. After 
each time step the prescribed pressure is updated from thin airfoil theory, thereby 
accountii^.for the displacement thickness interaction. Although this technique seems 
feasible, it needs further examination since Briley and McDonald obtained a regular solu- 
tion at a laminar separation point for a case with no interaction. The absence of the 
singularity in this case is probably due to numerical inaccuracy since a first-order 
scheme was used with a coarse grid. Hence, it is not clear in those cases in which 
interaction was included whether the solution at separation would be regular if a second- 
order scheme were used. 

In using the inverse boundary -layer procedures discussed above, it is necessary to 
incorporate a description of the inviscid flow to completely describe a viscous-inviscid 
intera'ctidrir ~ F6r~example , Kuhn and Nielsen (ref. 7) developed an iterative procedure in 
which Murman's inviscid transonic flow program (ref. 12) is solved iteratively with their 
inverse boundary-layer technique. Kuhn and Nielsen made calculations for the turbulent 
separated flow behind a bump placed on a wind-tunnel wall. The skin friction was updated 
for each new boundary-layer calculation based on the difference between the pressure 
computed from the boundary -layer solution and that obtained from the inviscid flow cal- 
culation for the displacement body. This procedure is not straightforward as it is not 
clear how to update the skin friction based on this pressure mismatch. A simpler pro- 
cedure, which was recommended by Kuhn and Nielsen and is the subject of the present 
paper, is to combine the displacement-thickness-prescribed, inverse boimdary- layer 
procedure with a suitable representation of the inviscid flow. In the present study of 
incompressible flow, the inviscid calculations are made by using an inverse Cauchy inte- 
gral from thin airfoil theory. 
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The Cauchy integral relating pressure and body slope is used in an inverse manner 
so that a new displacement body can be found from the pressure obtained in the inverse 
boundary- layer calculation. The displacement thickness is then updated by subtracting 
the prescribed body from the displacement body and the iteration is continued to con- 
vergence. This procedure is similar to that used by Jobe and Burggraf (ref. 13) and by 
Melnik. and Chow (ref. 14) in solving the asymptotic equations developed by Stewartson 
(ref.‘ 15) governing the flow at the trailing edge of a flat plate. An important feature of 
this inverse procedure is that the results require no smoothing, as has been typically 
found necessary in direct calculations of inviscid-.viscous interaction (refs. 16 and 17). 
This result is not surprising since the inverse procedure primarily uses numerical inte- 
gration, which is inherently a smoothii^ process, in contrast to the numerical differenti- 
ation used in the direct procedure. 

A second purpose of this paper is to present some calculations which have been 
made for turbulent flows involving separation and reattachment. These calculations, 
which do not include interaction, have been made by incorporating a two-layer eddy- 
viscosity model in the displacement-thickness-prescribed procedure used in the laminar 
analyses. The variable grid scheme analyzed by Blottner (ref. 18) is incorporated in 
these calculations to reduce the number of grid points across the boundary layer and 
simultaneously maintain second-order accuracy. 


SYMBOLS 


Cn>^n J 

CM 

Cf 

Cp 


coefficients in tridiagonal system of equations 

coefficients in Thomas algorithm 
skin-friction coefficient 
pressure coefficient 


L reference length 

£ coefficient in vorticity-transport equation (see eqs. (1) and (8)) 

m,n‘ indices for |- and 77-directions, respectively 


N computational coordinate 
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increment in computational coordinate 
free-stream Reynolds number, 
reattachment point 
separation point 
thickness (see fig. 3) 
free-stream velocity 
velocity component parallel to surface 

velocity component parallel to surface at edge of boundary layer 

Reynolds stress component 

coordinate along surface 

increment in coordinate along surface 

upstream interaction boundary 

downstream interaction boundary 

integration variable 

coordinate normal to surface 

coefficient of artificial time terna 

increment in displacement thickness 

boundary-layer thickness 

displacement thickness 



6 

eddy-viscosity coefficient 

r] 

transformed y- coordinate 

At7 

grid spacing in 77-direction 


molecular viscosity coefficient 

V 

kinematic coefficient of viscosity 

1 

transformed x-coordinate 

A| 

grid spacing in ^-direction 

P 

density 

riy 

stream function 




transformed stream function 

/ 

w 

vorticity 

aJ 


Subscripts: 


B 

body 

DB 

displacement body 

FP 

flat plate 

max 

maximum 

o 

denotes value at upstream boundary 

tr 

transition 



VISCOUS FLOW 


Governing Equations 

The governii^ equations for an incompressible boundary -layer flow for a prescribed 
displacement thickness can be given in terms of the followii^ vorticity-transport and 
stream-function equations: 



( 1 ) 




-♦2 _ 
= 6 (1 - T})ci> 


These equations are solved subject to the boundary conditions 


( 2 ) 


(3) 

- (4) 

where 6* is the displacement 
thiclmess which is defined in the usual manner. Equations (1) and (2) are nondimensional 
and the dependent variables which are barred have been scaled in the usual manner by 

\jR^ ^ which is appropriate for laminar boundary layers. The transformed stream 
function ^ is related to the usual stream function by xj/ = xf/ - ul (77 - 1) where u 
is the x-component of velocity. The vorticity is denoted by w and after it is obtained 
from equation (1), u is given by 

u(?,7?) = 6* yj (5) 


■ uUfi) = x}yU,0) = 0 

oJ(4,7j) and 0 as 77 — ®9 

The independent variables are | = x and 77 y/l* 


After equations (1), (2), and (5) are solved, the unknown edge velocity Ug is obtained 
from the x- momentum equation which is evaluated at the surface to give 



1 


6 * 


177=0 


( 6 ) 


This value should agree with the value obtained from equation (5) and thus .serves as a • 
check on the calculation. Equation (6) is integrated from the .upstream boundary to give, 
Ue(^) and the correspondii^ pressure coefficient is given by 
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(7) 


Cp = l-Ue2 

Further details of the preceding formulation for laminar flow are presented by Carter 
(ref. 6). 

For laminar flows the quantity £ appearing in equation (1) is set equal to unity, 
whereas for turbulent flow 


e = 1 + - (8) 

p 

where e/p is the ratio of the eddy viscosity to the molecular viscosity coefficient. The 
ed^-viscosity coefficient is used to relate the Reynolds stress to the velocity gradient 
in the usual manner 

= (9) 

In the present calculations a two-layer eddy-viscosity formulation has been used which is 
similar to that used by Harris (ref. 19), Cebeci and Keller (ref. 20), and others. In the 
inner region a combination of Prandtl's mixing lei^th model along with the Van Driest 
dampii^ factor is used and is given as follows in dimensional quantities: 


^ banner 


(0.4yD)2 

8U 

V 

8y 


( 10 ) 


where the damping factor D has been modified for separated flows and is given by 


D = 1 - exp 



1 

8u 


26 

V 

9y 

max 


( 11 ) 


In the outer region Clauser's velocity defect model is used along with the Klebanoff inter- 
mittency factor 


outer 


0.0168ue6* 

1 

V 

1 + 5.5(y/6)® 


(12) 


where 6 is the boundary-layer thickness defined as the point where u = 0.995ug. The 
boundary between the inner and outer regions is the point at which 
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(13) 


(f^lnner " (/"I 


outer 


Equations (10) to (13) relate the Reynolds stress to the mean flow and thereby complete the 
formulation for turbulent flows. This two-layer eddy-viscosity model has been widely 
used for attached flows with mild pressure gradients; its applicability to separated tur- 
bulent flows is unknown at the present time. Nonetheless, the purpose of the present 
paper is to develop a numerical scheme for separated turbulent flows; in the future, 
refinements of the turbulent empiricism will be made to assess the quantitative results 
in comparison with experiment. 

Numerical Procedure 

Figure 1 gives a schematic diagram of the computational schemes and boundary 
conditions used for the boundary-layer calculations. In an earlier paper Carter (ref. 6) 
presented a global iteration procedure in which the finite-difference scheme is switched 
in the reversed-flow region to properly account for the reversed-flow direction. The 
global iteration technique requires repeated streamwise iterations until convergence is 
obtained. More recently. Carter and Wornom (ref. 21) have shown that a separated bound- 
ary layer can be computed much more rapidly with the usual forward-marching procedure 
used for attached boundary layers provided that the tridiagonal equations are diagonally 
dominant and that the streamwise convection of vorticity is neglected in the reversed-flow 
region; that is, if u < 0, then set 

^ = 0 (14) 

This approximation is somewhat similar to that used by Reyhner and Fliigge-Lotz (ref. 22) 
for neglecting the streamwise convection of momentum in the reversed-flow region. Use 
of this approximation eliminates the well-knowji instability of marching in a direction 
opposite to that of the flow; in addition, the accuracy is essentially unaffected if the mag- 
nitude of the reversed-flow velocity is less than O.lUoo, as shown by Carter and Wornom 
(ref. 21). 

The streamwise gradients which are typically encountered as a boundary layer 
approaches separation are quite large, and thus, from numerical experimentation it has 
been found necessary to use a fully second-order -accurate scheme for u > 0. The 
C rank- Ni cols on scheme, including the reversed-flow approximation discussed previously, 
is used in the laminar calculations and has been previously discussed in references 6 
and 21. The turbulent forward- marching procedure, which is presented subsequently, is 
quite similar although the computational molecule shown in figure 1 is used since the 
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Crank- Nicolson scheme resulted in oscillatory solutions for some of the turbulent cal- 
culations. In this case the streamwise derivatives are represented by the second-order 
difference expression * , 

(15) 

The bars, which were used previously to denote Reynolds number scaling, are deleted for 
convenience in the presentation of the numerical procedure. The computational molecule 
for u < 0 in the global turbulent calculation is the same as that reported for laminar 
flow by Carter (ref. 3). A constant grid spacing is used in the stream direction and 
hence 4 = *« A|. In the normal direction a variable grid is used, second-order accuracy 
being maintained by using the difference expressions developed by Blottner (ref. 18) which 
are given as 


l^‘^m,n ■ ^‘^m- 


^ = ‘^m,n+l ~ ‘^m,n-l ^ q/^2) 

"^m,n Vl'Vl 


(16) 


0r?2 'Vl'Vl ^n+l-^n s ^n ’ Vl 

m,n I 


Blottner showed that a variable grid scheme is equivalent to a coordinate stretching - 
method if the coordinate 77 , can be related to a computational coordinate N in which 
the grid is evenly spaced. In the present calculations the grid is varied at a constant , 
rate; that is, = K which can be written in terms of a computational coordi- 

nate N as follows 


^n ^max 


Nn/AN 
K ^ - ] 

kVan . 1 



where = (n - 1)AN and 0 = N = 1. In the present calculations K, the ratio of adja- 
cent grid spacings, equals 1.09, ~ 21> ^*^*1 points are used across the bound- 

ary layer. This grid point distribution insures a minimum of 15 points in the viscous sub- 
layer. However, no numerical study was made to determine the optimum value of K or 
the minimum number of grid points. By using equations (15) to (17) the vorticity trans- 
port equation can be written in the following form where q denotes the column iteration 
level: 
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A9a.5r!. 1 + (Bn + , = dJ t o<nJ, „ 

n m,n-l \ “ ! m,n n m,n+l n m,n 

where ’ 

An = -(C^ + K£m n_i) 

Bn = (K + l)j(!m,n + 3C| 


Cn - Ctj - ^ni,n+l 


Dn = " ‘^m-2,n) 


a - 2jC^| + £ni,n+l + ^^m,n-l " (K + l)^ni,n 


In the reversed-flow region the streamwise convection of vorticity is neglected; as a 
result = 0. 

Repeated application of equation (19) from the wall to the outer boundary results in 
a tridiagonal system of linear equations for the vorticity. These equations are solved 
the Thomas algorithm which can be written as 

1 ^ (21) 

m,n n n m,n-l j ' / 

where . , . 


JJ. Pg - 

Bn + CnC„^i 


Bn + CnC^+1 


‘ ( 22 ) 


In equation (19) a timelike term a (“’mi - “’m,n) has been introduced to provide the 


unconditional diagonal dominance 
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• |Bn + a|^ |An| + lCn| . ^ ' ■ (23) 

and thereby prevent error growth in the back substitution procedure given in equation (21). 
Introduction of this term is a modification of the usual implicit technique used to solve 
the boundary-layer equations and is the subject of a recent paper by Carter and Wornom 
(ref. 21). 

The quantities bj^ and are computed, beginning at the outer boundary where 
the boundary condition o>(|,«>) = 0 is imposed and proceeding to the wall. Equation (21) 

is then used to obtain once the value at the wall is known. The wall vor- 

m,n m,l 

ticity is found by simultaneously solving for the stream fvmction from equation (2) across 
the boundary layer and imposing the surface boundary condition given in equation (3). 
Details of this procedure are the same as those used for laminar flow which are pre- 
sented by Carter (ref. 6). After the back substitution in equation (21) is completed, the 
coefticients in equation (19) are updated and the process continued until convergence is 
obtained. Convergence is assumed when the maximum change in all the dependent vari- 
ables between two successive column iterations is less than 10“ 3, In most of the calcu- 
lations it was necessary to use underrelaxation for the iterative column solution, as was 
discussed by Carter (ref. 6). The relaxation factor typically used in both the laminar and 
turbulent calculations was 0.6. ‘ 

INVISCID FLOW AND INTERACTION PROCEDURE 

In this section the inviscid flow, which is approximated by small-disturbance theory, 
and the interaction procedure shown in figure 2 are discussed. This procedure is applied 
to the laminar flow over the surface shown schematically in figure 3 in which flow sepa- 
ration and reattachment occur. The strong viscous-inviscid interaction region is 
assumed to be limited to the region shown in figure 3 in order to replace the infinite 
limits in the Cauchy integral, which is used to compute the inviscid flow, with finite val- 
ues. Thus, it is assumed that the region of strong interaction is located a large distance 
from the leading edge and the Blasius flat-plate displacement thickness at this location 
results in a negligible pressure gradient. These assumptions are discussed in further, 
detail later. 

The calculation is begun with an assumed displacement thickness distribution for 
Xq s X i Xj^, which is input to the boimdary-layer equations. The resulting solution yields 
the surface pressure Cp from equation (7) which results from the inviscid flow over 
the ..Replacement body. The displacement body coordinate, ygg is given by ^ 

yDB = yB + (24) 
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where yg is the prescribed body surface. The displacement thickness can be written 
as 


5*(x) 


1.7208Vx 
+ A(x) 




(25) 


where the first part is the Blasius flat-plate solution and the second part A(x) is the 
result of the nonzero pressure gradient. Since the inviscid flow is linear, the pressure 
on the displacement body can be written as 


" S,B ^ S,FP 


(26) 


where Cp 3 is the pressure on the prescribed surface when no viscous effects are pres- 
ent and is found from the direct Cauchy integral 



where the u^al small- disturbance approximations have been made. It is necessary, of 
course, to compute Cp B only once. In equation (26) Cp^A is the pressure coefficient 
due to A, the deviation of the displacement thickness from that generated by a flat plate. 
Also, in equation (26) Cp^FP denotes the pressure coefficient induced by aTlat-plate 
displacement thickness, which in the present study will be approximated as 

CpFP = 0 (28), 

This approximation, which is discussed by Van. Dyke (ref. 23), is the result of second- , 
order boundary -layer theory for the flow over a semi -infinite flat plate. Furthermore, 
it should be noted that the pressure gradient induced by a flat-plate displacement thick- 
ness increases as the leading edge is approached and thus, if Cp fp * then there is no 
logical point at which a downstream, relatively local interaction calculation can be initi- 
ated other than at Xq = 0. . 

Since 'the boundary-layer solution is computed inversely, the iteration procedure is 
simplified by also solving the inviscid flow with an inverse technique. The inverse, 
Cauchy integral relating the pressure on the displacement body to the rate of growth of 
the displacement body is given by 
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^5B(x)=l.r 

dx 2ff ‘^-oo X - X’ 


(29) 


Similarly, the inverse relation for the prescribed surface is given by 

dx 2it^-co X - X' 


(30) 


If equation (30) is subtracted from equation (29), then from the previous discussion it fol- 
lows that 


dx 2it J-oo X - x' 


(31) 


The end points on this integral pose a problem since the boundary-layer solution 
which gives 1® computed only in the range Xq ^ x S Xj. In the present study the 

lower limit on the integral in equation (31) has been replaced with Xq and thus, the inter- 
action is assumed negligible upstream of Xq| that is, = 0 for x^Xq. Downstream 

of the interaction at x = xj, it is expected that Cp^^ will be small but not zero, partic- 
ularly in the early stages of iteration. It is noted that since the prescribed surface 
returns to a flat plate as x — «>, then Cp ^ = 0 is the required downstream 

boundary condition. This boundary condition is imposed in the inviscid calculation since 
the boundary-layer problem is parabolic and is solved independently of a downstream 
boundary condition. To avoid the discontinuity which would be encountered by setting 
Cp,A = 0 at X = xi, the following extrapolation is used for x = xj: 





(32) 


where Cp^A denotes the extrapolated value of Cp^^ and the coefficients a^, a.2, and 
a 3 are computed by matching equation (32) with Cp^A obtained in the boundary-layer 
solution at X = xj. Numerical tests on the approximate treatment of the limits in equa- 
tion (31) will be discussed later. 
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On the basis of the preceding discussion, equation (31) is rewritten as 



Cp,A(x')dx' ^ r°° Cp-A(x')dx’ 

X - x' X - X' 



The first integral is evaluated numerically with a third-order-accurate quadrature formula 
in which the singularity at x' = x is isolated in the same manner as was done .by Jobe and 
Burggraf (ref. 13). The second integration is performed analytically with equation (32). 

In solving equation (33) dA/dx is obtained at the midpoints between the grid nodes of the 
boundary -layer solution. Thus, A at the boundary-layer points x = | = m A^ is then 
obtained to second-order accuracy from , . , 

^m = Am_l+^ Ax + 0(Ax2) (34) 

beginning at x = x^ where A = Aq which is found by solving the boundary-layer equa- 
tions in Gortler variables from the leading edge with Cp^B prescribed. Note that Aq 
must be small since the interaction is assumed negligible upstream of x = Xq. As shown 
in figure 2 the new displacement thickness is computed frpm equation (25), and the result- 
ing vailue is multiplied by ^oo,L to give . 6* = ^oo,L®* order to conform to the 

usual boundary-layer scaling. At this point in the iteration cycle a check on convergence ■ 
is made which is defined in the present study as, 

max < 10"^ (35) 

m 

where i denotes the iteration cycle. Generally, it is found that when equation (35) is 
satisfied, the corresponding maximum change in Cp is about 10"®. Underrelaxation.was 
used in the present calculations at the indicated points in figure 2. A relaxation factor 
of 0.2 was used in the calculations; several attempts were made to increase this value, . 
but these calculations diverged. , ■ 

RESULTS AND DISCUSSION 

Turbulent Boundary- Layer Calculations 

Calculations with the inverse boundary-layer procedure are first discussed for the 
turbulent flow over a flat plate in which the input displacement thickness is, computed from 
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a direct solution of the boundary-layer equations in Gortler variables. A comparison of 
the skin-friction distribution obtained from’ the direct and' inverse calculations is shown in 
figure 4 along with the experimental measurements of Wieghardt and Tillman (flow num- 
ber 1400 in the 1968 Stanford Conference,, ref. 24). The excellent agreement between the 
direct and inverse boundary -layer calculations is expected since the same eddy-viscosity 
model is used in both of these calculations. In the inverse calculations it is found that 
with Ax = 0.025 the deduced edge velocity deviates from unity by 5 percent, whereas with 
Ax = 0.0125 this error is less than 2 percent. These boundary-layer solutions are also 
shown in figure 4 to be in good agreement with the experimental data. The direct solution 
starts at the leading ec^e with a laminar boundary layer and transition is assumed to occur 
at X = 0.001m. The inverse calculation begins at x = 0.087m, the direct solution at that 
station being used for the upstream boundary condition. Better agreement with the data 
could be obtained either by using a more detailed modeling of the transition region or by 
using the experimental data for the starting conditions, as was done by most of the inves- 
tigators in the Stanford Conference. 

Further calculations were made by usii^ the displacement thickness distributions 
shown in figure 5 as input conditions to the present inverse boundary -layer procedure. 

The resulting solutions for the skin friction and ec^e velocity are shown in figures 6 to 8. 
Figure 6 shows the large separated region computed for laminar flow with the 6 * dis- 
tribution designated as case 1 in figure 5. Note that the edge velocity, and hence the pres- 
sure, shows a plateau region between separation and reattachment which is characteristic 
of separated flows.. This laminar calculation is the same as that designated as case B in 
reference 21 where additional details are presented. The turbulent boundary -layer solu- 
tion corresponding to this same 6* distribution is shown in figure 7, where it is seen that 
separation did not occur despite the 30-percent decrease in the edge velocity; The quali- 
tative trend here is correct since it is well known that a turbulent boundary layer requires 
a larger pressure rise than a laminar boundary layer before separation occurs. As a 
check on the solution shown in figure 7, the deduced velocity Ug was, used as an input to 
the boundary -layer equations expressed in Gortler variables. The resulting skin-friction 
distribution is seen in figure 7 to be in good agreement Avith that obtained in the inverse 
solution. There is some difference in the two solutions near the point of majdmum 
boundary-layer thickness and is probably due to a lack of resolution across the boundary 
layer in the direct solution. 

A more severe case was computed by using the displacement thickness distribution 
designated as case 2 in figxire 5 as the prescribed condition. The resulting skin-friction 
distribution is shown in figure 8 and it indicates that separation and reattachment occur 
for this case. It is seen that there is no difference between the solution obtained with the 
approximate forward- marching procedure and that found by the global iteration technique, 
in which the finite-difference scheme is switched in the reversed-flow region to properly 
account for the reversed-flow direction. In this case the magnitude of the reversed-flow 
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velocity is only about 0.02Uoo, which at least for laminar flows is quite accurately com- 
puted by the approximate forward- marching procedure, as shown in reference 21. 

An attempt was made to compute the same case as a direct calculation with Ue, 
prescribed to see whether or not a singularity exists at the separation point as in laminar 
flow. The preliminary results indicate that the direct solution has a much steeper skin- 
friction gradient near separation, although further grid refinements as well as analysis • 
are needed to determine the precise behavior in this region. Furthermore, the signifi- 
cance of this study is unclear since the solution behavior near separation will depend on 
the turbulent formulation, which is quite approximate in the present paper. 


Laminar Viscous-Inviscid Interaction 

The prescribed surface, for which interaction calculations are presented, is given 
by 


yg = t sech 4(x - 2.5) ^ (36) 

where t = -0.03 in figure 9 and t = -0.015 in figure 12. For these results Xq = 1.0, 

Xj = 4.0, and Ax = 0.025 which results in 121 grid points in the x-direction. The 
Reynolds number based on free-stream conditions and the distance from the leading edge' 
to the assumed start of the interaction is = 8 x 10^. In the boundary -layer calcula- 

tion 87 points were used across the viscous layer; thus, a total of 10 527 grid points were 
used which were found to require approximately 25 sec on the CDC 6600 computer to com- 
plete one iteration cycle shown in figure 2. The results shown in figure 9 were obtained 
with a relatively crude initial guess on 6* as seen in figure 10 and were found to require 
64 iterations to meet the convergence criterion given in equation (35). No attempt was 
made to optimize the convergence rate of these calculations. 

In the lower part of figure 9 the deduced displacement body is shown in comparison 
with the prescribed body. The points S and R refer to the separation and reattach- 
ment points, respectively, and are connected by the dividing streamline which separates 
the inner recirculating flow from the outer, forward flow. The displacement thickness 
distribution is better seen in figure 10 where the initial and final distributions are com- 
pared with the Blasius flat-plate solution. A comparison of the pressure distributions on 
the displacement and prescribed body shapes is shown in the upper part of figure 9. The 
difference in these two curves is which is about half of the uninteracted pressure 

level at the bottom of the trough; this shows the large influence the boundary layer exhibits 
in this flow field. It is observed that approaches zero as both the upstream and 

downstream boundaries are approached, which indicates that the interaction is adequately 
contained in this region. Additional calculations were made to insure that the results are 
independent of Xq, xj, and the extrapolation given in equation (32), provided that these 
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boundaries are placed sufficiently far from the bottom of the trough. As a further numer- 
ical test the grid spacing in the x-direction was reduced from 0.025 to 0.0125; the results 
of these two calculations differ only slightly. From previous studies (ref, 6) the grid 
spacing across the boimdary layer is considered sufficiently small and thus, the results 
showii in figures 9 and 10 are considered an accurate solution of the governii^ partial 
differential equations. 

The skin-friction distribution for this case is shown in figure 11 along with the 
Blasius flat-plate distribution and that obtained by a direct calculation of the boundary- 
layer equations with the inviscid pressure distribution prescribed. This latter calcula- 
tion demonstrates the usual singularity at separation where the slope of the skin-friction 
curve is infinite, whereas the solution obtained including interaction has a large but finite 
gradient at separation. The column iterative procedure used in the direct calculation no 
longer converges downstream of the separation point. As expected, the point of separa- 
tion is predicted further downstream when the effects of interaction are included since the 
boundary layer reduces the adverse pressure gradient and thus delays separation. It is 
observed that the gradient is even larger near reattachment and is followed by an over- 
shoot of the flat-plate result which is characteristic of the usual neck region downstream 
of a separated flow. Comparison of figures 10 and 11 shows that the maximum in skin 
friction corresponds to the minimum in displacement thickness in the neck region. 

The results of an additional calculation in which t = -0.015 in equation (36) are 
presented in figures 12 and 13. Comparison of these results with those discussed pre- 
viously for t = -0.03 show that the shallower trough results in a smaller but not negli- 
gible viscous-inviscid interaction. In figure 13 it is seen that for this case the inclusion 
of the viscous-inviscid interaction relieves the adverse pressure gradient such that the 
flow remains attached, despite the prediction of separation from a first-order boundary- 
layer calculation using the inviscid pressure distribution which is given in figure 12. 

CONCLUDING REMARKS 

In the present paper a technique is demonstrated for solving laminar and turbulent 
separated boundary layers. The turbulent separation results appear qualitatively correct; 
however, the quantitative accuracy of this solution technique must be evaluated by making 
comparisons with experimental data. It is anticipated that modifications of the eddy- 
viscosity model used in the present study will be required. 

A new viscous-inviscid interaction procedure for separated flows is discussed and 
several calculations using this technique are presented. This iterative procedure requires 
none of the smoothing techniques usually required in numerically matching a boundary- 
layer and inviscid flow solution. Further studies are needed to optimize the efficiency of 
this interaction procedure. 
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Figure 1.- Boundaxy-layer computational schemes and boundary conditions. 



Figure 2.- Inverse viscous-inviscid interaction procedure. 
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Figure 3.- Schematic diagram of interaction region. 
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Figure 4.- Flat-plate skin-friction comparisons. 



Figure 5.- Prescribed displacement thickness. 








Figure 8.- Skin friction for case 2 - turbulent separation and reattachment. 
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